Spin-dependent inelastic collisions in spin-2 Bose-Einstein condensates 
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We studied spin-dependent two-body inelastic collisions in _F = 2 *^Rb Bose-Einstein condensates 
both experimentally and theoretically. The ^'^Rb condensates were confined in an optical trap and 
selectively prepared in various spin states in the F — 2 manifold at a magnetic field of 3.0 G. 
The measured atom loss rates depend on the spin states of colliding atoms. We measured two 
fundamental loss coefficients for two-body inelastic collisions with total spins of and 2. The loss 
coefficients determine the loss rates of all the spin pairs. The experimental results for mixtures of all 
spin combinations are in good agreement with numerical solutions of the Gross-Pitaevskii equations 
that include the effect of a magnetic field gradient. 

PACS numbers: 03.75.Mn, 34.50.-s, 03.75.Kk 



I. INTRODUCTION 

The field of cold collisions has attracted extensive in- 
terest and has grown explosively since the early days of 
atom cooling and trapping Jj] . The development of novel 
techniques for cooling and manipulating atoms have led 
to a deeper understanding of physics of collisions: e.g., 
evaporative cooling enabled Bose-Einstein condensates 
(BECs) to be realized and opened the field of ultracold 
collisions, while the evaporative cooling process itself re- 
lies on the nature of collisions. The mean-field theory 
of BEC, which successfully describes the properties of 
BEC, depends on the elastic s-wave scatteringlength to 
characterize the atomic interaction energy Inelas- 
tic collisions also play important roles, e.g., in cooling 
and retaining an atomic or molecular cloud in a trap 

In an optical trap, the spin degrees of freedom of 
atoms are liberated enabling a rich variety of spinor BECs 
physics to be studied. Spin-2 BECs have attracted much 
interest in recent years and there have been several ex- 
perimental studies of spin-2 systems, including investi- 
gations of their inagnetic phases [1, 0, H, [§| , of multiply 
charged vortices and the phase separation between 
spin-2 and spin- I BECs It is, however, difircult 

to observe stationary states of these phenomena because 
the ®^Rb F — 2 BEC has a negative scattering length, 
and in "^Li, ^^Na, ^^K, ^"^Rb the = 2 states are in the 
upper hypcrfine manifold. Ultracold gases in the upper 
hyperfine states decay to the lower hyperfine state. *^Rb 
atoms are highly suitable for investigating the properties 
of a spin-2 BEC because they have much lower inelas- 
tic collision rates than other species: e.g., ^^Rb has a 
loss rate of the order of lO^^^'cm'^/s, which is much lower 
than that of ^SNa 0. 

In this paper, we systematically investigate ultracold 
two-body inelastic collisions between Zeeman states of 
F = 2 atoms both experimentally and theoretically. The 



experiment was conducted by creating an = 2 spinor 
^^Rb BEC in an optical trap and populating desired 
Zeeman sublevels by radio-frequency transitions. We 
measured the atom loss rates for various initial popu- 
lations. Spin-exchange collisions were negligible at a bias 
magnetic field of 3.0 G. By analogy with the scattering 
length in elastic collisions, two-body inelastic collisions 
are described by two parameters, bo and 62, which cor- 
respond to channels with the total spins of and 2, re- 
spectively. We experimentally determine these two pa- 
rameters from the loss rates of single-component BECs 
of \F = 2,mF = 1) and \F = 2,mF = 0). We also 
measured the atom loss rates for two-component BECs 
of all possible sets of magnetic sublevels and we com- 
pared the results with our theoretical model. We calcu- 
lated the time evolution of the number of atoms in two- 
component BECs using the values of 60 and 62 obtained 
by two methods: the single-mode approximation (SMA) 
and the Gross-Pitaevskii (GP) equations that include the 
effect of a magnetic field gradient. The results obtained 
with the latter method agreed well with the experimental 
results. 



II. THEORY OF TWO-BODY INELASTIC 
COLLISIONS 

We first consider a system of spin-2 *^Rb atoms with 
no atom losses. The Hamiltonian for this system is given 

by 



^. n 



where V'm and Vm are the field operator and an external 
potential for the hyperfine state \F = 2, m), and M is the 
mass of the ®^Rb atoms. The interaction Hamiltonian 
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ffint has the form UA 
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with ayr being the s-wave scattering length with the col- 
hding channel of total spin and Ajrj^ is the annihila- 
tion operator of two atoms with total spin |JF, ^A) defined 
as 
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FIG. 1: (color online) Inelastic collision rates /s'j"'''" ^ for all 
combinations of m and m' in a spin-2 BEG. A rectangle that 
connects two ellipses corresponds to a collision between two 
different components, whereas a rectangle linked to one ellipse 
corresponds to a collision between the same component. 



with C^:^, = {T,M\2,m;2,m') being the Clebsch- 
Gordan coefficient. 

We take into account only the two-body inelastic loss, 
since the two-body inelastic collision is the dominant 
atom loss process for the upper hyperfine state of *^Rb 
atoms. We neglect the three-body losses and density- 
independent losses, such as background-gas scattering 
and photon scattering. The two-body inelastic collision 
changes the hyperfine spin of either or both of the col- 
liding atoms from = 2 to = 1. Consequently, the 
colliding atoms acquire a kinetic energy corresponding to 
the hyperfine splitting, allowing them to escape from the 
trap. 

We write the master equation for the time evolution of 
the system with two-body inelastic loss as 



dp 
dt 



(5) 



loss 



where p represents the density operator for the atoms in 
the F — 2 manifold. The second term on the right-hand 
side of Eq. (O describes the atom loss, which must be 
rotationally invariant. If we assume that the atoms after 
the two-body inelastic collision immediately escape from 
the system and do not interact with the other atoms, the 
atomic-loss part of the master equation can be written 

as fii 



dp 
dt 
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where the constant hj: characterizes the loss rate for col- 
liding channel with total spin T . We note that Eq. ([6]) is 
rotationally invariant. If the dipolar decay is negligible, 
the total spin T before and after the collision must be the 
same, and therefore, inelastic decay through the = 4 
channel is prohibited in Eq. ([S]): this is experimentally 
confirmed in the following section. 



For example, let us consider a mixture of |2,0) and 
|2, 1) atoms. Using Eq. ([S]), we obtain the two-body in- 
elastic loss of the 1 2, 0) component as 



= Tr 
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where Tr represents the trace. In the mean- field approx- 
imation, Eq. ([7]) is simplified as 



where K2^'^^ — bo/5 + 2b2/7, K'^'^' = 62/7, and is the 
density of component m. In general, two-body inelastic 
decay in a two-component mixture has the form 

——^^K^ ■ >n„^{r)-K'2 ^n„(r)n„, (r), (9) 

and for a single-component system, 
drimir) 
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dt 
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The coefficients K. 
Fig.ffl 



(m.m') 



for all m and m' are shown in 



From Eq. ([5]), the equation of motion for is given 



by 

^(V'™(r)) = ^([V;™(r),i?])+Tr 



' loss 
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Using Eq. (I6|), the second term on the right-hand side of 
Eq. (fTT|) can be rewritten as 
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We find that Eq. (fT2|) has the same form as ([-ffint, "0™]) 
with gjr being replaced by hjr/2. In the mean-field the- 
ory, we thus expect that the two-body loss can be in- 
corporated into the GP equation by replacing gjr with 
gj= — itibjr/2. Thus, the GP equation for a spin-2 BEG 
with two-body loss is obtained as 



ih 



dt ~ V 2M 
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where 
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50 = 30 



.92 = 92 



54=34. (14) 



If we assume that the atom density rim(r, t) of a single- 
component BEG always follows the Thomas-Fermi pro- 
file. 
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where am is the scattering length for component m and 
nit) is determined from 
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Nm{t) = J nmir,t)dr, 
we can integrate Eq. ([T0| , giving [3, [§] 



dt 
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Here the coefficient is given by 
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with Lu being the average trap frequency. The solution of 
Eq. (fT7|) is obtained as 
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5/2 ■ 



(19) 



For a two-component system, we assume that both com- 
ponents m and m' have identical density profiles that are 
given by 



nm{r,t) 
nm'ir,t) 



-^p„(t)[M(i)-nr)], (20) 
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where d is the arithmetic average of the three scattering 
lengths (between m-m, m'-m' and m-m'), Pm(m')(t) = 
Njn{m'){t)/N{t), and fi(t) is determined from 

Nit)= [ [nmir,t)+nm'ir,t)]dr. (22) 



We also assume that the coherent spin dynamics is neg- 
ligible. Such an SMA is accurate if the three scattering 
lengths are the same. Integrating Eq. ([9|), we obtain 



dNmjt) 
dt 
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where 7 has the same form as Eq. (|18p in which a„j is 
replaced by d. 



III. EXPERIMENTAL SETUP 

The experimental apparatus and procedure to create 
^^Rb condensates are almost the same as in our previ- 
ous study [1, H]. A precooled thermal cloud prepared 
by the first magneto-optical trap (MOT) was transferred 
to the second MOT by irradiating the first MOT with a 
weak near-resonant cw laser beam focused on the center 
of the MOT. Over 10^ atoms in the second MOT were 
optically pumped into the |2,-|-2) state and recaptured 
in an loffe-Pritchard (clover-leaf) magnetic trap. A BEG 
containing 10^ atoms was created by evaporative cooling 
with frequency sweeping an rf field for 18 s. 

The BEG was loaded into a crossed far-off-resonance 
optical trap (FORT) at a wavelength of 850 nm. One 
of the laser beams propagated along the long axis of the 
cigar-shaped BEG and had a power of 7 mW and a 1 /e^ 
radius at the focus of 31 /im. The other laser beam had 
a power of 11 mW and a beam waist of 97 /im and was 
focused on and crossed the first beam. The potential 
depth of the crossed FORT was estimated to be about 
1 /iK. Both lasers had power fluctuations of less than 
1%. The radial (axial) t rap frequency measured from 
the parametric resonance [1^ was 2tt x 162 Hz (27r x 21 
Hz) . The average trap frequency determined by measur- 
ing the parametric resonance was 2it x 82 Hz, which was 
consistent with that determined by measuring the release 
energy [l7| . 

The condensed atoms in the |2,+2) state in the mag- 
netic trap were transferred to the crossed FORT. The 
magnetic trap was quickly turned off and the quantiza- 
tion axis was then nonadiabatically inverted. The con- 
densed atoms were kept in the trap for 200 ms. In the 
way, condensates containing 3.2 x 10^ atoms in the |2, —2) 
state were prepared in the crossed FORT. Some of the 
condensed atoms in the |2, — 2) state can be transferred to 
desired spin states |2,tof) by the Landau-Zener transi- 
tion with an external magnetic fleld of 20.5 G [l^ . In this 
magnetic field, the condensed atoms can be selectively 
transferred to the desired target states by quadratic Zee- 
man shifts. The relative populations of \2,mp) states 
were controlled by the rf field strength and the sweep 
rate. The observed fluctuation in the number of atoms 
was estimated to be 10%, which was due to variations in 
the initial number of BEG atoms and in the rf-transfer 
rate between spin states. After state preparation, the 
20.5-G fleld was turned off and an accurately controlled 
external magnetic fleld of 3.0 G was immediately applied. 
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Trap time (ms) 

FIG. 2: (color online). Time dependence of the number of 
atoms in a single component BEG initially prepared in |2,0) 
(filled circles), |2, — 1) (open squares), and |2, — 2) (open tri- 
angles) states. The inset shows the decay of |2, —2) atoms for 
a longer period. The experimental data are fitted by Eq. (|17p 
(dashed lines). The theoretical curves are calculated from the 
GP equation (solid lines). 



The crossed FORT was turned off at the end of the time 
evolution period, and absorption imaging was performed 
after a free expansion time of 15 ms to measure the 
population distribution of each spin component by the 
Stern-Gerlach method [l^ ■ A magnetic field fluctuation 
of ~10 mG was measured by observing the transfer rate 
for magnetic dipole transitions and the residual gradient 
magnetic field was estimated to be ^30 mG/cm Q. 



IV. RESULTS AND DISCUSSION 

Figure [5] plots the time dependence of the number 
of atoms of a single-component BEG initially prepared 
in the |2,0) (filled circles), |2,— 1) (open squares), and 
|2,— 2) states (open triangles). Each point represents an 
average typically over four samples. The observed con- 
densate fraction was greater than 90% and the tempera- 
ture of the thermal cloud was below 100 nK. In the mea- 
surement, finite temperature effects in the condensates, 
as seen in Ref. [l^, were not observed. The relative pop- 
ulation of the initially prepared single component was 
observed to be larger than 0.95 through the time evolu- 
tion since a magnetic field strength of 3.0 G suppressed 
the spin-exchange collisions and only the initially pre- 
pared single component was observed through the time 
evolution [^. The quadratic Zeeman energies between 
|2,0) and |2,±1) and between |2,0) and |2,±2) were 31 
nK and 124 nK, respectively. The number of atoms in 
|2,0) and |2, —1) states decreased rapidly with time. By 
contrast, the number of atoms in the |2,— 2) state re- 
mained almost unchanged within 100 ms, showing good 

( — 2 —2) 

agreement with the collision rate if 2 ' =0. The 
value of the one-body loss coefficient Ki is experimen- 
tally estimated to be Ki = 0.4 s~^ from the result for 
the |2,— 2) state (see the inset of Fig. [5]), and the three- 



body rate constant K3 was deduced by Soding et al. Q 
to be = 1.8 X 10~^^ cm^s~^. In our experimental 
condition, the contributions of the Ki and A'3 terms to 
the atom loss are negligible compared with that of the 
K2 term. Therefore, the decays in the number of atoms 
for the |2,0) and |2, -1) atoms are well fitted by Eq. fTT)) 
based on SMA calculation. The inelastic collision rates 
are thus obtained as icf '"^ = (8.9 ± 0.9) x lO-^^^cm^/s 
and K^2^^'^^^ = (10.4±1.0) x lO-i^^cm^/s for = and 
mp = —1, respectively. From these values of K'f''^^ — 

60/5-I- 262/7 and ATj"^'"^^ — 3^2/7, we deduce the inelas- 
tic collision coefficients 62 = (24.3 ± 2.4) x 10~^''cm^/s 
and 60 = (9.9±5.5) x 10^^^ cm'^/s with scattering lengths 
of 94.57aB for mp — Q and 95. 6803 for mp = — 1 with 
a_B being the Bohr radius, respectively [§]. We also cal- 
culated the time evolution using the GP equations with 
the above values to compare the results obtained using 
the GP equations with those obtained using the SMA. 
The difference between the results obtained by GP and 
SMA calculations is negligibly small as shown in Fig. [21 
We next consider the loss of a two-component mixture 
to check the consistency of Fig. [T] combined with the mea- 
sured values of and 62- We observed atom losses for 
all pairs of spin states with initially equal populations as 
shown in Fig. [3l The atom loss rates depend strongly on 
the spin states of colliding atoms. We first discuss the 
atom loss rate for the (2, — 1) + (2, —2) pair. As Fig. El^a) 
shows, the number of |2, —2) atoms remained almost un- 
changed even in the presence of |2,— 1) atoms. This is 

f — 1 —2) ( — ^ —2) 

consistent with the prediction K2 ' ~ = 

as shown in Fig.[T] On the other hand, for (2, 0)-|-(2, —2), 
(2, +1) + (2, -2), and (2, -1-2) + (2, -2) pairs, the |2, -2) 
component decreased with time as shown in Figs. [3] (b), 
(c), and (d), respectively. This decrease can be inter- 
preted by two-body inelastic collisions of |2,— 2) with 
|2,0), |2,-|-1), and |2,-f2) components because of the 
large values of if^+^'"^' = 662/7, AT^ '"^^ = 462/7, and 
is:^+2'-2) ^ 45^/7 260/5 from Fig. [H The dotted and 
dashed lines are the SMA calculations from Eq. us- 
ing the measured values of 60 and 62, an average trap 
frequency of 27r x 82 Hz, and an average scattering length 
between spin states. The calculations are in qualitative 
agreement with the experimental results for all pairs. In 
particular, for pairs without the |2,— 2) component (i.e. 
the (2,0) + (2,-1) and (2,-^1) + (2,-1) pairs shown 
in Figs. [3] (e) and (f), respectively) the SMA calcula- 
tions with the Thomas-Fermi approximation are in good 
agreement with the experimental results indicating the 
prediction of Fig[TJ 

However, with regard to the |2,— 2) component, the 
SMA calculations of the (2, 0)-h (2, -2), (2,-M)-h(2,-2), 
and (2, -1-2) + (2, —2) pairs are slightly smaller than the 
experimental results as shown in Figs. [3] (b)-(d). The 
decay curves for (2, -hi) + (2,-1) and (2,-^2) + (2,-2) 
pairs are theoretically predicted to be identical since 
is:^~^'~^^-|-X^+^'"^^ = 462/7-^260/5 is equal to A:^+^'"^\ 
This is also understood from the fact that |2, -|-2)-t-|2, —2) 
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FIG. 3: (color online). Time dependence of numbers of condensed atoms for various initially populated spin states (a) 
(2, ~1) + (2, -2), (b) (2, 0) + (2, -2), (c) (2, +1) + (2, -2), (d) (2, +2) + (2, -2), (e) (2, 0) + (2, -1), and (f) (2, +1) + (2, -1) 
pairs. The theoretical curves are calculated from Eq. (|23|l with the SMA (dotted and dashed lines) or the GP equations that 
include the effect of a magnetic field gradient (solid and dot-dashed lines). 



is transformed to |2, +1) + |2, —1) by spin rotation. The 
SMA calculation results for (2, +1) + (2, —1) pair are in 
good agreement with the experimental results. However, 
the calculation for the (2, +2) + (2, —2) pair is smaller 
than the experimental results for trap times longer than 
30 ms. 

Figure m shows relative center-of-mass position in axial 
direction (z-axis) of BECs between different spin states. 
The center-of-mass position of mp = m atoms for a trap 
time t, Zmit), was measured after Stern-Gerlach separa- 
tion. The separation between the center-of-mass position 
of mi? = m and mp = rn' atoms due to the Stern-Gerlach 
magnetic field is 2;„i(0) — Zm' (0) because two components 
overlapped at t = 0. The relative center-of-mass po- 
sition for a trap time t shown in Fig. 3] is given by 
[zm{t) - Zm'{t)] - [zm{0) " Zm'{0)]. The magnetic field 
gradient could be the reason for the difference in the rela- 
tive center-of-mass positions between different spin states 
since the low-field- and high-field-seeking states move in 
opposite directions. The relative center-of-mass position 
for the (2, -1-2) -I- (2, —2) pair shifted to positive values 
with time while that for the (2, -|-1) -I- (2, —1) pair leveled 
off almost zero. The separation could reduce the inelastic 
collision loss of the (2, -1-2) -I- (2, —2) pair. 

The separation may also be caused by mean-field in- 



teractions of BECs. The miscibility of multicomponent 
BECs is determined by s-wave scattering lengths between 
identical and different spin states. The interaction energy 
for a mixture of spin m and m' states is given by 



E-mt = I dr 



5n 



where gn 



(24) 

4:Trhamm'/M with amm' being the scat- 
tering length between m and m' states, amm' can be 
expressed in terms of a^r by comparing Eq. (j24p and the 
mean-field expectation value of Eq. Q for a mixture of 
TO and m' states pT| . Two-component condensates is 
immiscible when gmm9m'm' 



< 9 r, 



Therefore, based 



on the measured scattering lengths given in Ref. [8| the 
(2, 0) -I- (2, -1) and (2, -1) + (2, -2) pairs are expected 
to have a phase separation while the other pairs are ex- 
pected to be miscible. 

Figure[5]shows optical densities in the axial direction (z 
axis) of all spin-state pairs at a trap time of 100 ms, which 
were averaged in the direction of gravity. These density 
distributions depend on the spin states of the pairs. Bi- 
nary BECs of (2,-1) + (2, -2) and (2, 0) + (2, -1) pairs 
exhibit phase separation and domain formation as shown 
in Figs. [51^ a) and (e), due to the immiscibility condition. 
The (2, 0) -I- (2, -2) and (2, +1) + (2, -2) pairs exhibit 
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FIG. 4: (color online). Relative center-of-mass positions of 
(a) (2, 0) + (2, -1) pair (open circles) and (2, 0) + (2, -2) pair 
(closed circles), (b) (2, +l) + (2, —2) pair (closed triangles) and 
(2, -l) + (2, -2) pair (open triangles), (c) (2, +l) + (2, -1) pair 
(closed triangles) and (2, +2) + (2,-2) pair (open triangles). 



some inhomogeneity between different spin states. The 
(2, +2) + (2, —2) pair exhibit phase separation despite 
it being predicted miscible on the basis of Eq. and 
some other pairs have partial phase separation or homo- 
geneous distribution in Fig. [5] This behavior could be 
explained by the effects of the magnetic field gradient. 

The above analysis demonstrates that it is necessary 
to consider both the effect of the magnetic field gradient 
on motions in the FORT and interactions between BECs 
in the theoretical analysis. We numerically solved the 
GP equation with V„i = Vtrap + B' ^B'm'Z/2, where 
the field gradient B' in the axial direction was assumed 
to be 30 mG/cm, the radial (axial) trap frequency was 
27r X 162 Hz (27r x 21 Hz), and the spin-dependent scat- 
tering lengths were derived from Ref. The results 
are shown by the solid and dot-dashed lines in Fig. [31 
These theoretical results are good quantitative agreement 
with the experimental results. The differences between 
the theoretical and experimental results are significantly 
reduced for the (2,-1) + (2,-2), (2,0) + (2,-2), and 
(2, +2) + (2, -2) pairs as shown Figs, [^a), (b), and (d). 
In particular, the calculation for the (2, +2) + (2, —2) 
pair can describe reduced atom losses after 50 ms and 
a kink around 30 ms as shown in Fig. El^d). However, 
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FIG. 5: (color online). Axial profiles of averaged density 
distributions for various spin states at an evolution time of 
100 ms. Initially populated to (a) (2,-1) + (2,-2), (b) 
(2,0) + (2,-2), (c) (2,+l) + (2,-2), (d) (2, +2) + (2,-2), 
(e) (2, 0) -I- (2, -1), and (f) (2, +1) + (2, -1) states with equal 
population. 



the measured displacement of the relative center-of-mass 
position of the (2, -t-2) -I- (2, —2) pair cannot be explained 
by magnetic field gradient alone. We do not have an 
adequate theory for describing it quantitatively. In the 
(2,-1-1) -I- (2,-2) pair, the theoretical resuh for |2,-|-1) 
atoms cannot completely explain the experimental re- 
sults, while that for |2, —2) atoms is in better agreement 
with the experimental results. This may be caused by 
the inhomogeneous density distribution. The lower atom 
losses than the SMA calculation is due to the separation 
of the two components and the suppression of two-body 
collision between them. 

It has been theoretically predicted Q and experimen- 
tally measured [1] that the magnetic phase of a spin-2 
^^Rb BEG is in close vicinity to the boundary between 
the cyclic and antiferromagnetic phases. Saito and Ueda 
[2l| proposed a method to distinguish these two phases. 
If the magnetic phase is cyclic, a (2, -1-2) -f (2, —2) mix- 
ture produces (2, 0) -I- (2, 0) pairs, and observation of the 
|2, 0) atoms in the (2, -1-2) -I- (2, —2) mixture may be a sig- 
nature of the cyclic phase. However, the spin-dependent 
two-body losses are not taken into account in Ref. [2l[. 
To determine the magnetic phase for the upper hyper- 
fine state, it is necessary to take into account the spin- 
dependent inelastic collision rates @. Even if the |2,0) 
atoms are produced by the cyclic interaction, the inelastic 
collision channels of the (2, 0)-F(2, -^2) and (2, 0)-H(2, -2) 
states may eliminate the |2, 0) atoms, thereby erasing the 
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evidence of the cyclic phase. The obtained coUision rates 
would enable us to quantitatively predict the production 
rate of |2,0) atoms for the cyclic phase. Another diffi- 
culty for determining the magnetic phase is the inhomo- 
geneity of the system. In Ref. ^] , it is assumed that 
all condensates overlap homogeneously over the whole 
evolution time. However, in the measurement of the 
(2, +2) + (2, —2) pair, two components were separated 
in time evolution, as shown in Fig. [5Ud). The magnetic 
phase in the ground state could be determined by spin 
population measurement when condensate separation is 
prevented by controlling the magnetic field gradient or 
applying an additional external field such as optical lat- 
tices. 



V. CONCLUSIONS 

We have studied ultracold two-body inelastic collisions 
between Zeeman states of F = 2 atoms both experi- 
mentally and theoretically. We measured spin-dependent 
atom loss rates for all possible combinations of spin states 
by creating F = 2 spinor ®^Rb BECs in an optical trap 
and populating the desired Zeeman sublevels by radio- 
frequency transitions. Atom loss rates were measured at 
a bias magnetic field of 3.0 G at which spin-exchange 
collisions are negligible. We show that the two-body in- 
elastic collision rates between all Zeeman sublevels can be 
expressed in terms of two inelastic collision constants, bo 
and 62, which respectively characterize collision channels 



with total spin (60) and total spin 2 (62). We have deter- 
mined these values to be 62 = (24.3 ± 2.4) x 10~^'*cm^/s, 
and bo — (9.9 ± 5.5) x lO"^** cm^/s from measured de- 
cay rates of single-component BECs (i.e. the |2,0) and 
|2,— 1) states). Observed time evolutions of the atom 
number for two component BECs were consistent with 
SMA calculations using the values of bo and 62. The 
discrepancy between the SMA calculations and the ex- 
perimental results can be explained by spatial separation 
between two components. The calculated results from 
the CP equations with a magnetic field gradient were in 
good agreement with the experimental results. 

A detailed understanding of the relative-population de- 
pendence and the spin-state dependence of inelastic col- 
lisions are key issues in the future study of spinor BECs, 
such as the determination of the magnetic ground state 
of spin-2 ®^Rb BEC and observation of novel quantum 
vortices. 
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